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Abstract 

We generalize the concept of conserving, <I>-derivable, approximations to rel- 
ativistic field theories. Treating the interaction field as a dynamical degree 
of freedom, we derive the thermodynamic potential in terms of fully dressed 
propagators, an approach which allows us to resolve the entropy of a relativis- 
tic plasma into contributions from its interacting elementary excitations. We 
illustrate the derivation for a hot relativistic system governed by electromag- 
netic interactions. 



*It is a pleasure to dedicate this paper to Leo Kadanoff on his sixtieth birthday. Immediately after 
Leo and I received our Ph.D.'s from Harvard in 1960, we began working together at the Institute 
for Theoretical Physics in Copenhagen (now the Niels Bohr Institute) on the problem of how 
to construct self-consistent approximations to two-particle propagators that preserved the basic 
conservation laws. Our solution to the problem was written up in the paper Conservation laws 
and correlation functional, and the concept generalized to deriving selfjconsistent approximations 
from a functional, <£, of the single particle Green's function in Ref.cl. I hope that Leo enjoys 
revisiting these early ideas, which may prove useful in understanding modern problems of the 
thermodynamics and transport properties of systems with long-ranged gauge fields. - Gordon 
Baym 
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I. INTRODUCTION 



The motivations for studying relativistic plasmas are numerous. Nuclear matter at very 
high densities is expected, from the asymptotic freedom of QCD, to be in the form of a 
plasma of deconfmed quarks and gluons ||. Such plasmas were present in the early universe 
and may exist in the cores of neutron stars; current experiments aim to produce and study 
them in collisions of ultrarelativistic nuclei M. Electromagnetic plasmas are also of interest 
as Abelian models of quark-gluon plasmas. In addition, problems of relativistic plasmas are 
closely related to current issues in condensed matter theory. In gauge fields models of high 
T c superconductors and of the fractional quantum Hall effect f|, the magnetic component of 
the interaction plays a fundamental role; this component is suppressed in normal metals by a 
factor (vp/c) 2 (where vp is the Fermi velocity), but is important in both strongly correlated 
systems and relativistic gases. 

The thermodynamics and quasiparticle modes of relativistic plasmas have been the sub- 
ject of much earlier work. In order to understand the equation of state of interacting rela- 
tivistic plasmas, the thermodynamic potential has been calculated up to the first six orders 
of perturbation theory in the coupling constant . Studies of the plasma microscopic 
properties have also revealed the existence of well-defined quasiparticle excitations and col- 
lective modes []5|JT0[] . In particular, the fermionic spectrum has certain remarkable features 
not encountered in non-relativistic plasmas: the spectrum has a gap at zero momentum, 



and splits into two branches at small momenta [JTO 



Many physical quantities of relativistic plasmas are infrared divergent when evaluated 
in perturbation theory, as a consequence of the lack of static screening of magnetic inter- 
actions. Taking into account dynamical screening of such interactions at long wavelengths 
(the anomalous skin effect [JIT]), equivalent to including effects of Landau damping in the 
photon polarization operator, eliminates the divergences in transport coefficients ]12| and in 
the rate of energy loss of fast particles Polarization effects also modify the interaction 
between low energy quasiparticle modes. The "Hard Thermal Loop" expansion scheme, 
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proposed by Braaten and Pisarski |14| , handles the effects of the medium diagrammatically 



by introducing vertex and self-energy corrections at small four-momentum. This description 



can also be cast in terms of kinetic equations 15 



Even after taking into account corrections from medium effects, difficulties still remain. 
Due to a lack of static screening in the magnetic component of gauge interactions, the 
fermion quasiparticle damping rate at finite temperatures is still divergent in perturbation 
theory. Blaizot and Iancu have shown that |T(| in QED at finite temperatures, including 



multiple scattering a la Bloch-Nordsieck in the vacuum leads to well-defined, divergence- 
free, quasiparticle modes. This structure of the fermion propagator also appears in gauge 
field models of the fractional quantum Hall effect ||17|| . However, it is not clear whether 



the infrared structure in the fermion spectrum actually affects physical observables such as 
transport coefficients and the specific heat |I8|| . Even quantities that are well-controlled 
order by order in perturbation theory can cause difficulties. The (asymptotic) expansion 
of the thermodynamical potential converges very slowly, which has lead to suggestions to 
reorganize the perturbation expansion in terms of dressed fermion states rather than those 
of a free interacting gas JT9[ . 

Our aim here is to give a general framework for analyzing the effects of different char- 
acteristics of the fermion spectrum, the presence of a gap and the infrared structure, on 
thermodynamic quantities. Our analysis is based on ^-derivable conserving approxima- 
tions, first introduced P,|20[| in the context of quantum transport theories. Among various 



applications, their use in the study of liquid 3 He has allowed one to interpret the T 3 log T 
term in the low temperature specific heat in terms of repeated scattering of particle-hole 
pairs pT] , p2| . The technique has been extended to Bose condensed systems in Refs. p3 



A generalization of the conserving approximation techniques to relativistic plasmas at zero 
temperature provided a controlled expansion of the ground state energy of electromagnetic 
and quark plasmas up to order g 4 |7j . In this latter work, the free energy is given as a func- 
tional of the fully dressed fermion and boson propagators as well as fully dressed vertices. 
In our present analysis, we keep track of both fermion and boson elementary excitations 



and treat both the matter and interaction fields as dynamical quantities. This allows us to 
decompose the entropy in the form 




where the entropy densities of free fermion and boson gases, Of = — /log/ — (1 — /) log(l — /) 
and at = —n\ogn+ (1 +n) log(l +n) are here weighted by the spectral densities A s and B s 
of the interacting system. 

In this paper, we develop ^-derivable approximations for a hot relativistic QED plasma, a 
gas of electrons and positrons in a thermal bath of photons at temperature T, where T ^> m, 
although we expect the approach to be valid also for non-Abelian theories including QCD. 
For simplicity, we set the electron mass, m, to zero. We then illustrate the technique by 
computing the entropy within the one-loop approximation, and comment briefly on the 
general structure of the entropy spectral densities A s and B s . 



II. THE THERMODYNAMIC POTENTIAL 

We first review ^-derivable approximations in non- relativistic field theories 0. The 
thermodynamical potential of a non-relativistic Fermi system can be written as a functional 
of the fully dressed fermion propagator G as, 

/3Q = +trSG-tr log(-G). (2) 

Here (3 = 1/T and the symbol tr denotes a trace over spins and coordinates. In four- 
dimensional notation, 

trX^Y, [~ t0 dlX(l,l + ), (3) 

er,tr' 

where 1 is shorthand for the coordinates (ti,rx), while 1 + indicates that the time variable 
in the second argument of X is t\ + i0 + . The functional 3>[G] is defined as the sum of 
all skeleton graphs contributing to the diagrammatic expansion of the potential Q, where 
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all lines are fully dressed propagators, G, instead of bare ones, Go- In contrast with the 
familiar diagrammatic expansion of a scattering diagram in the vacuum, the expansion for 
$[G] contains a symmetry factor 1/n for each diagram with n fermion lines. Upon variation 
of G, the functional $[G] satisfies the relation 

8$[G} = tr£5G, (4) 

where £ is the electron self-energy. Pictorially, this relation means that removing one of 
the n fermion lines in a given ^-diagram produces a self-energy diagram with n—1 fermion 
lines. Since there are n ways of removing a line, no symmetry factor appears on the right 
side of Eq. (H). The self-energy £ is here a functional of the fully dressed propagator G, 
and it satisfies Dyson's equation G -1 = Gq~ 1 — £, where the bare propagator Go is given by 
Go 1 = 5 a ^i{uj n — p 2 /2m) for a non-relativistic system of fermions. 

The principle behind conserving approximations is the following: one selects a particular 
subset of diagrams for the functional, $ a [G], from which one deduces an approximate self- 
energy functional S a = 5§ a [G}/5G. Dyson's equation provides then a self-consistent relation 
for the approximate propagator G a . As shown in ||, the approximation for G leads to current 
densities and an energy-momentum tensor that obey the continuity equations expressing 
charge, particle number, energy and momentum conservation. The basic ingredient of the 
approximation scheme of Eq. (fj) which enforces conservation laws is the stationarity property 
of Q, 5Q = 0, upon a variation of G that keeps Go constant, as can easily be checked from 
Eqs. @ and (§. 

We now generalize this approach to relativistic field theories, illustrating the method for 
electromagnetism, where the Lagrangian is 

C = $(i$-e4)iJ>-±F'»'F tu ,. (5) 

In QED, the free energy f3VL is given in terms of the fully dressed electron and photon 
propagators G and D, by 

(3tt = $[G, D] — Tr EG + Tr log(- 7 °G) + ^TrlED - ^Tr log(-L>) = W[G,D}. (6) 
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Here, the functional $[G, D] is the sum of all the skeleton diagrams of the thermodynamic 
potential, expressed in terms of fully dressed G and D instead of the bare electron and 
photon propagators, Go and Dq. Under a simultaneous variation of the propagators G and 
D, 

5$[G, D] = Tr ESG - ^Tr USD, (7) 

as one easily sees by removing an electron line or a photon line in a given diagram contribut- 
ing to $. The electron self-energy T,[G,D] and the polarization operator U[G, D] are here 
functionals of G and D and satisfy Dyson's equations 

G- 1 = Go 1 - S, (8) 
D- 1 = Dq 1 - n. (9) 

To prove that the functional D] in Eq. (|6]) is identical to the thermodynamic 

potential (30. of an electromagnetic plasma, we scale the coupling by a parameter A in the 
interaction Lagrangian, 

Ant = (10) 

We work in Coulomb gauge. In the non-interacting limit, A — > 0, both E and II vanish, and 
we recover the result for a gas of free electrons and photons 

(3Q = W[G , D ] = Tr log(- 7 oG ) - Tr log(-Avr), (11) 
where only transverse polarizations (subscript T) appear in the photon contribution [51 . 



To prove that (30, = W[G, D] it suffices to show that 

dQ x _dW 

~d\ ~ "9A~ (12) 

for arbitrary A. 

Consider first the quantity !3dVl\/d\. From the definition of the thermodynamical po- 
tential, 
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Q x = -Tlogtr [exp (-/3 (if A - //Q))] (13) 

(where Q is the charge operator, the Hamiltonian is H\ = H — J d 3 x C- m t, and tr denotes a 
sum over all the quantum states of the system,) we find 

^ = T {H * ) = Al* x{Clab) - ( 14 ) 

Here, (X) denotes the thermal average tr [exp(— (3{H\ — nQ)) X]/tr [exp(— (3{H\ — fiQ))]. 
We evaluate next dW[G\, D\]/d\. A given diagram contributing to $>[G\,D\] depends on 
A implicitly through the electron and photon propagators, G\ and D\, but also explicitly 
through the coupling Ae, Eq. (|ToD . However, the terms arising from the dependence on Ae 
through G\ and D\ cancel in the partial derivative, since VK[Cr, D] is stationary under a 
variation of G and D. From the definition of D], the variation of W is 

SW[G, D] = TrS5G - ^TtU5D - Tr5(SG) + ^-Tr5(IW) + Tr 5{\og{-j G)) 

-iTr5(log(-D)). (15) 

Using 51og(— 70G) = —GbG~ x and Dyson's equation, SG^ 1 = — <5£, as well as similar 
relations for the photon propagators, we find that 

8W[G, D] = 0, (16) 

under a variation that does not affect Go? Dq, or the coupling constant, Ae. Therefore, in 
the derivative of W with respect to A only the explicit dependence in the coupling constant 
contributes, and we find 



dW <9$ 

-[G\,D\] = ——[G\,D\] 



(17) 

G.D 



d\ 1 ' AJ d\ 

where the derivative on the right side is taken at fixed G and D. 

Next, to compute the right side of Eq. fllTf) , we use the invariance properties of $. As 
each vertex of a given ^-diagram is connected to two electron lines and one photon line, the 
functional $ remains constant under the scaling transformations 
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$[G A , D\] = ®[s- f G x , s b D x ; sX], f = 1 + 6/2, 
$[G A , D A ] = $[s-^G? A , s l D x] A], / = 6/2, 



(18) 
(19) 



where the third arguments of $ on the right sides are the scaling factors of the coupling 
constant. Taking the derivative of Eq. (|I8"D with respect to s, one finds 



A 



<9$ 



d(s\) 

which, for s = 1, gives 



G,D 



-/^- 1 Tr^G A + 6^Tr§D A = 0, 



A 



G.D 



f Tr S A G A + - Tr U X D X = / Tr S A G A + (/-!) Tr n A Z) A . 



Similarly, the derivative of Eq. fll9|) with respect to s gives 



TrE A G A + Trn A D A = 0. 



Therefore, Eqs. (H), (0) and (0) give 



1 

^rT-[G A ,D A ] = — TrS A G A . 
cm A 



(20) 



(21) 



(22) 



(23) 



To complete the proof, we relate <9Tr£ A £ A /<9A to (H int ), starting with the equation of 
motion for the time-ordered fermion propagator G A (1, 1') = — i(Tijj(l)ijj(l')) , which, from 
the Lagrangian, Eq. (|), is 

0iG x (l, 1') + zAe(T4(l)^(l)^(l')) = 5(1 - 1'). (24) 



Comparing to Dyson's equation in real space, 

J dl G^{1, 1)G X {1, 1') -Jdl £ A (1, T)G A (T, 10 = 5(1 - 10 



(25) 



we identify 



TrS A G A = Jdl JdlX x (l,l)G x (l,l + )=iXeJ dl $(1)^1(1)^(1)) 
= (3Xe [ dr 1 (^(r lj 0)4(r 1 ,0)V(ri,0)) = / 9(/r ilIt >. 



Therefore, 



dW 1 /3 dfi 

/3 aA-=A TrSAGA= A ( ^ t)=/5 9A' 



(26) 



(27) 



which completes the proof. 



III. THE ENTROPY 



We now turn to the derivation of the entropy of the plasma in terms of the fully dressed 
propagators G and D, which is given by the derivative 



on 

b ~ df 



= -^(TW[G,D]) 



(28) 



at constant volume V and chemical potential fi. We work in frequency-momentum space, 
where the electron and boson propagators have the spectral representations 

G(w n ,p) = / , (29) 

D T {u n , q) = / — i q - , (30) 

DLM = \ + fP^^, (31) 



q 2 J 2-k uj n — uj, 



where the subscripts L and T denote longitudinal and transverse polarizations. The Mat- 
subara frequencies are uo n = m(2n + 1)T + \x for electrons and uo n = 2mnT for photons. The 
functional D] depends therefore on the temperature through the spectral functions A 



and B, through the complex frequencies u n , as well as through factors T at the vertices |25 
Because D] is stationary under a simultaneous variation of the propagators G and D, 

Eq. fljjp , it is also stationary under variations of A and B that ignore the other temperature 
dependences. Thus, the contribution to the entropy coming from the temperature deriva- 
tives of the spectral densities alone vanishes. In the following, we therefore evaluate every 
derivative with respect to the temperature at constant A and B. 

We start from the expression of the thermodynamical potential j3VL = W[G, D], Eq. (^j). 
Evaluating the frequency sum by standard contour integration techniques [p5|1 , we find 

du, 



{l = T<5>iG,D] + Y, I -^/K)trlm \ZG + log(- lo G~ 

p J 7T 

+ E -f^r n M E 9i im [n, a + M-A" 1 )] . ( 32 ) 

q J A7t l=L,T 

where f(u> p ) = (exp(P(uj p — //)) + and n{uj q ) = (exp(/3uj q ) — 1) _1 are the Fermi and Bose 
occupation factors, tr is a trace over spinor indices, the degeneracy factors gi are = 1 for 



longitudinal modes and qt = 2 for transverse ones. Here, Im[F], means lm[F(u + ?0 + )]. 
The integration over u q is a principal value integral, since in deforming the contour integral 
onto the real axis, one needs to go around the (Bose) frequency u> n= o = 0. Differentiating 
Eq. ( |32| ) with respect to T at constant A and B, we decompose the entropy into a sum of 
three terms, 



on 

b ~ df 



= S f + S b + S', (33) 

fi,V,A,B 



where 



*/ = -E/^|tr (imS ReG + Im log^C" 1 )) , (34) 
Sb = -Y.9ij^^ (indl, ReA + Im log(-Df')) , (35) 



S' 



d{T<$>) 



dT 



E/^ftr(ImGReE) 



A,B 

-E^/^^I-ARen,. (36) 

In this decomposition, the terms Sf and Sb are the contributions from the electron and 
photon elementary modes, while 5" is a correction term arising from the interactions. 

We illustrate the derivation in the one-loop approximation, which corresponds to taking 
the diagram $ depicted in Fig. 1, and the self-energy diagrams of Fig. 2. Applying the usual 
Feynman rules, and carrying the sum over frequencies, we find 

167r V p^,l J UJ p + UJ q - UJ p > 

(37) 

= Tl^Eft / dujqdup'Siicjp, ,p')Bi{u q ,q) — — - ^ + - - + ^ + S HF , (38) 

where / = f(u) p ), f = f{u)p>) and n = n(u q ), z is complex, p' = p + q, and the matrix 
elements P\ and S\ are given by 
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P L (u p ,p;uj p ,,p') = tr 7°A(u;p,p)7 A(^p/,p / ) 



P T (u p ,p;uj p ,p') = -(S i:j - quq,)tr \j t A(u p ,p)y'A(u p r,i/ ) 
S L (u p ,,i/)='fA(u p ,,i/)'f, 
S T {u p >,p') = -(S i:j - qi^fA^pf,^)^. 



(40) 
(41) 
(42) 
(43) 



The terms $hf and S H f in Eqs. fl37|) and (|38"D arise from the static term 1/q 2 in the photon 
propagator of Eq. (|3~T|) , 



HF 



J HF 



8tt 2 1/ 



^ / diOpduJpi tr 7°A(u;p,p)7 A(u;p/,p / ) 



27rV ^7 

p 



^ / £2uy7°A(uy,p')7° 4 



(44) 
(45) 



In the one-loop approximation 5" in fact vanishes. (A similar result was observed in the 
SPA approximation in liquid 3 He by Riedel [22|.) To see how the various terms in Eq. fl3"S] ) 
cancel, we first note that the same combination of matrix elements appears in each term of 
S'. Since lm[G{uj p + i0+)} = -A(u p ,p)/2 and Im[A(wp + iO + )] = -B l {u q ,q)/2, 



S' 



8n 2 V 



p,p' 



^ / dujpdupf tr [7°A(^p,p)7 v4(uy,p' ) 



<ShF 



+ 



16vr 3 V 



9ij~ dujpdup'dujg Pi(uJ p ,p; uj p/ ,p)Bi(uj g , q) 



S 



UJ V + UJ a — UJ V > 



(46) 

where the S and iShf denote combinations of statistical factors and their T derivatives. 
Using Eqs. (|3~7| ) through (fi5|)) symmetrizing the contributions of ReEnF and ReS by the 
transformation (u> p ,p) <->(ay, p'), (<^ g ,q) <->(— co> g , — q), and using the fact that the spectral 
functions E>i are odd functions of their arguments, we have 

Shf = ~d {ff} /dT + fdf/dT + fdf/dT = 0, (47) 
S = -d {n/(l - /') - (1 + n)(l - /)/'} /AT + (/(l - /') - /'(l - /)) dn/dT 
+ (n(l - /') + (1 + n)/') + (n/ + (1 + n)(l - /)) d/'/dT = 0. (48) 

Thus the correction term S' vanishes in the one-loop approximation. 

This cancelation takes place only in the lowest order diagram for $. A general analysis 
of $ diagrams in the context of non-relativistic normal Fermi liquids |21| shows that the 
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contributions to 5" come only from those graphs that have at least two vanishing energy 
denominators. These terms correspond graphically to the set of diagrams that can be cut in 
three and only three different pieces when one removes a set of fermion lines. Generalizing 
this analysis to the diagrams representing the functional D], we see that there is only 
one possible cut in the diagram of Fig. 1 and it generates the vanishing denominator 
uo q + uj p — ujpi. The simplest diagram that contributes to S' is shown in Fig. 3, where we 
have drawn a set of cuts that give two vanishing energy denominators. 



IV. INTERPRETATION OF THE ENTROPY FORMULA 

A. Exchange and correlation entropy 

When the two electron lines in Fig. 1 are replaced by bare propagators and the photon 
line is replaced by a dressed propagator in the RPA approximation, i.e., with the lowest order 
bubble diagram insertions, we recover the set of diagrams that contribute to the exchange 
and correlation terms considered by Akhiezer and Peletminskii |J. It is instructive to see 
how the expression of the entropy that we have derived, Eqs. (Bl), fl3"5|) and (P§), generates 



these exchange and correlation terms correctly when expanded in the first few powers of the 
coupling constant e. We start with the exchange entropy ||, 

S, = ~{TV n«D ) = £ / („ Im n« D a ) , (49) 

where E^ 2 ) is the photon self-energy of Fig. 2 with two bare electron propagators. [In this 
section, we do not write the sum over polarization states explicitely.] Expanding Eqs. (p4|) 
and (|35|) to first order in e 2 , we find 

^=E/^^trImG ReS[D ]+E/^§Im^oRen( 2 ). (50) 



Then expanding Eq. ( pop for S' = to order e 2 , we recognize that the right side of Eq. (|50] 
is 5(2) = -d(T$W)/dT, which from the identity ^ = (-l/2)TrIl( 2 ) D becomes S® = 
(l/2)<9(TTrn( 2 ) D )/dT = S s . 
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The sum of the exchange and the correlation terms that contribute to the entropy is 
1 d 



S c + S s 



(tTt log(-l + D U {2) 

1** (to 

^ T 2n \dT 



log(-l + D U 



(2)> 









— n Im 




) 


dT 





(51) 



where D is now the propagator in the RPA approximation, D^ 1 = Dq 1 — iT 2 ) . Using 
arguments similar to those we used to derive the identity S' = 0, we can decompose the last 
term on the right side into the following terms: 



x - f dLJ q 

^ J 2tt 



n Im 



dT 



? /^| ImGoRe2|B| 

doj n On 



^^Imn (2) ReL>. 
2tt dT 



(52) 



The two terms on the right side correspond respectively to Sf expanded to linear order in 
and to the first term in Sb- Adding the first term on the right side of Eq. (|5lD, we have 
therefore shown that our expression for S = Sf + gives the correct correlation term in the 
entropy when one uses bare electrons and the RPA photon propagator in Eqs. (]3~4])- (|36D . 



B. Decomposition in elementary excitation modes 

The formulae we obtained for the entropy terms Sf, Eq. fl34|), and Sb, Eq. (p5|), allow us 
to perform a spectral analysis by decomposing Sf and Sb into integrals over the elementary 
excitations of the matter field and the electromagnetic field. We closely follow the derivation 
that was carried out by Carneiro and Pethick for liquid 3 He [ET]. The key is to transform 
the temperature derivatives of the statistical factors / and n into the following expressions 

df _ dof 



dT' duj p [UJph 



dn 
df 



(53) 
(54) 



where again <jf = —/log / — (1 — /) log(l — /) and a n = —nlogn + (1 + n) log(l + n) 
are the entropy contributions from an electron mode of energy u p and from a photon mode 
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of frequency u> q . Integrating Eqs. ([34]) and (|35|) by parts, we see that Sf and Sb obey the 
spectral representations 

s f = J2j ^VfMMupiP)* ( 55 ) 
S b = J2{ ^-Vb(u q )B s (uj q ,q), (56) 

q JO Z7T 

where A s and £> s are defined by 

A s (u p ,p) = tr A {ReGT + 21m log(- 7o G)} , (57) 
5 s (u; 9 , g ) =^^^-{ReA^ + 2Imlog(-A)}, (58) 

while r(dj p ,p) = — 2 Im T l (u p +iO + , p) and Li(u q ,q) = —2ImU[(u q +iO + , q) are the imaginary 
parts of the electron and photon self-energies. In deriving Eq. (p6[) , we have used the fact 
that the boson propagator D is an even function of its arguments to reduce the domain of 



integration in Eq. (|56|) to positive frequencies u g . 

The sum over spinor indices (the trace tr ) in Eq. ( p|) can easily be decomposed into two 
contributions, one from electron states with a chirality equal to their helicity (subscript +) 
and one from states with opposite helicities and chiralities (subscript — ), by writing 

q = 7o ~ 7 • P G+ + 7o + 7 • P q 



which gives 



trReGr = 2(ReG+£ + + ReGL£_), (60) 



where the factor 2 is from the spin sum. Hence, A s = 2 ^± A s ±. 

We now turn to the entropy spectral functions A s and B s . In the one-loop approximation 
of Fig. 1, the self-energies £ and II of Eqs. fl5B] ) and ( PP| ) depend on the fully dressed spectral 
functions A and B. These functions depend in turn on the self-energies £ and IT through 
Dyson's equations, Eq. (P). The problem of determining A and B, and the functions A s 
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and B s , is therefore a self-consistent one. However, we can make general statements about 
their structure on the basis of the following arguments. If the system develops well-defined 
excitation modes, we expect the functions A and B to consist of narrow peaks at the locations 
of the quasiparticles and collective modes, as well as wide bands of continuum states. The 
continuum states that contribute to the spectral function A are composed of particle-photon 
states, while those contributing to B are particle-hole and particle-antiparticle states. In 
lowest order of perturbation theory, the continua extend over frequency ranges — p < uo p < p 
and — q < u> q < q, respectively PJlQ|. We expect that in the present self-consistent problem, 
interactions only slightly modify these frequency ranges. The quasiparticle modes are at 
frequencies well outside the spectrum of continuum states, \u p \ > p and u q > q. 

The structure of the entropy spectral functions A s and B s is qualitatively similar to that 
of A and B. The photon spectral function B s , as we see from its definition, Eq. (|58|), has a 
support over the frequency range of the continuum states contributing to B. It also contains 
a contribution from the oscillation modes of the electromagnetic field, i.e., the longitudinal 
and transverse plasmon modes, with frequencies 0Ji{q) and u>T{q) arid momentum q, where 
Wl,t(q) > <?• As we show, the spectral function B s takes a simple form in the vicinity of 
these frequencies 

{ZiUf/2 . 
((uj q - ui{q)y + (ZiLi/2yy 

where Z\ = dReDf 1 / du q , Z L ~ uj L /{2q 2 ) and Z T ~ \/{2ut)- Evaluating the logarithmic 
term in Eq. (|58|), we have for each polarization state I, 

2Imlog(-A) = 2nQ(-ReDf x ) - 2 arctan — ^— T , (62) 

2ReZ); 

where the values of the arc tangent are in the range [—tt/2, tt/2]. The function Li/(2uj q ) is 
the interaction rate of an excitation of frequency u q . In the vicinity of the modes, u q ~ 
1jJ l,t{q)i Li/2uj q is perturbatively smaller than the modes frequencies uj LjT (q) f2H[| , and varies 
slowly around u L>T {q). Thus, with ReA = ReDf 1 / ((ReDf 1 ) 2 + {Li/2) 2 ), and ReDf 1 ~ 
{uj q — u>i(q))/Zi, we can evaluate the right side of Eq. Q5%p keeping Li constant, and derive 
Eq. ©. 
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The spectral functions B s i are therefore sharply peaked around the mode frequency u>i. 
Comparing to the photon spectral densities Bi, which have a Lorentzian form close to the 
poles, Bi ~ (ZiLi)/((ujq — uji) 2 + (ZiLi/2) 2 ), we see that B s j have a stronger peak and smaller 
wings. The electron spectral function, A s = 2 J2± A s ±, has a structure similar to that of B s : 
continuum states contribute at small energies |iG> p |^p, while A s ± exhibit a sharp variation 
close to the quasiparticle energies uo±, 

A s± (uj p ,p) ~ ^rtjP^ ( 63 ) 

{{up - u Q ±y + {z±r±/2yy 

where Z± is the quasiparticle mode residue. In deriving Eq. (|63D, we have neglected the 
variation of the (one- loop order) interaction rate r±, which for energies u v in the vicinity of 
uo± is ~ 0{e 2 T\og{qD/\u p — ct>o±|))> where qn ~ eT is the Debye momentum fll6"| . 

Both the electron and the photon excitation spectra exhibit collective modes at small 
momenta. However, the low energy modes contribute to the total entropy in very different 
orders in the coupling constant. The contribution from long wavelength modes of the elec- 
tromagnetic field has been calculated by Akhiezer and Peletmniskii |J ; the result is that the 
sum of plasmon oscillations and continuum modes contributes to the correlation term in the 
entropy, of order e 3 T 3 . We can easily understand this order of magnitude by noting that 
the contribution to the entropy of all states of energies of order of the plasmon frequency 
ul,t(0) ~ eT is ~ J2q ^ Qc-> where the momentum cutoff is q c ~ eT. The small energy 
modes in the electron spectrum give a smaller contribution, of order e 5 T 3 . To see this, we 
expand to lowest order the difference of Sf from its expression for a non-interacting gas, and 
concentrate on the phase space of small energy modes, uj p ~ eT, 

ASj ~ 2^E / ^^o ± K,p)ReSL 2) K,p). (64) 

For small energies, u p ~ eT, df/dT is ~ u;p/(4T 2 ), and A ±(u p ,p)~Re11±\up,p) is 
~ ±27cm 2 5{uj p z fp)/p, where irtf is the gap at zero momentum. Thus, ASf ~ rn 2 /T 2 J2 P ~ 
e 2 pl ~ e 5 T 3 , for a cutoff momentum p c ~ eT. The reason why AS/ is higher by two orders 
in the coupling constant than the contribution to the total entropy from long wavelength 
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modes of the electromagnetic field is the following. An expansion in power series of e of the 
logarithm term in B s , Eq. fl58|), leads to infrared divergences term by term 0; hence, to 
evaluate the correlation term in the entropy, one must include the logarithm as it stands. 
However, an expansion of the logarithm term in A s , Eq. (|57|), is divergence- free and results 



in a prefactor m 2 JT 2 ~ e 2 . 



Although the one-loop approximation is mathematically tractable, it has important lim- 
itations. To take into account correctly the low energy part of the electron spectrum, which 
enters only in order e 5 T 3 , one needs to consider ^-diagrams of higher orders, such as the one 
depicted in Fig. 3, contributing at least in order e 4 . Also, the diagram of Fig. 1 does not 
always provide the correct width for the entropy spectral functions evaluated close to the 
quasiparticle modes. For instance, the lifetime of a photon mode of energy ~ T is limited 
by Compton scattering and inverse pair annihilation. The self-consistent approximation of 
Fig. 1 includes only the direct terms of these processes, as can be seen by considering an 
electron self-energy insertion in one of the fermion lines of the bubble diagram shown in 
Fig. 2 and then taking the imaginary part of the corresponding photon self-energy. The 
cross terms arise from the imaginary part of the photon self-energy obtained by removing 
one photon line in Fig. 3. Finally, in calculating the widths of low energy modes, one also 
needs to include vertex corrections, as described by the "Hard Thermal Loop" perturbation 
scheme, see Ref. fL4| . 

We conclude by stressing the fact that our analysis of the entropy provides a framework 
for studying effects of the infrared structure in the fermion propagator on thermodynamical 
quantities. Kim et al. have proposed, in the context of the fractional quantum Hall ef- 
fect [0 , that one should be careful in calculating thermodynamical quantities by summing 
over fermion degrees of freedom, since by doing so, one may encounter infrared divergent 
terms. Kim et al. attribute this problem to the fact that the composite fermion propa- 
gator is not a gauge invariant quantity and suggest that the thermodynamical quantities 
should be calculated instead by summing over boson degrees of freedom only. In the present 
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analysis, the dangerous term is the first component of Sf, ~ J2 P J (du>p/2ii)df/dT Re GT, 
with T oc e 2 Tlog(g£ ) /|co'p — p\), and it does not lead to a divergence, as the principal part 
f-dx log | x | jx vanishes. This term is actually finite and is part of the exchange entropy 
term, of order e 2 . On the other hand, for energies u p in the vicinity of p, both the spec- 
tral densities A ~ T/((u p - pf + (r/2) 2 ) and A s ~ T 3 /2/((cu p - pf + (r/2) 2 ) 2 vanish as 
~ 1/r ~ log -1 \ojp— p\, as lu p ~ p. This logarithmic behavior is symptomatic of a breakdown 
of perturbation theory. In a future publication we will examine corrections in the free energy 
similar to those considered in the calculation of the fermion lifetime in a hot plasma [16|. 
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FIGURE CAPTIONS 

Figure 1: $ in the one-loop approximation. 

Figure 2: a) One-loop electron self-energy, b) One-loop photon self-energy. 

Figure 3: ^-diagram to explicit order e 4 . The dashed lines illustrate the cuts that give rise 
to a non- vanishing term in S'. 
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